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Abstract. The identification of very high energy 
photons is of great importance for the understanding 
of the origin of extreme energy cosmic rays (EECR). 
Several can be the sources of high energy photons at 
Earth. A guaranteed component is the flux of high 
energy photons expected as a consequence of the 
interaction of cosmic rays with the cosmic photon 
background. Another contribution may be expected 
as by-product at the acceleration sites of protons 
and nuclei, although such flux should be strongly 
suppressed for distant sources. On the other hand, 
top-down scenarios involving the decay of super 
heavy relic particles or topological defects, even 
if not currently favored, have as a characteristic 
signature an increasingly dominant flux of photons 
at the highest energies. In this work we study the 
statistical separation between hadron and photon 
showers at energies where both, LPM effect and 
magnetospheric interactions are important for the 
development of the cascades. We consider a detector 
with the same orbital characteristics as JEM-EUSO, 
but disregard trigger and reconstruction efficiencies, 
in order to define the maximum ideal discrimination 
power attainable. 

Keywords: extreme-energy cosmic rays; photon 
fraction 

I. Introduction 

The cosmic ray energy spectrum must have at least 
a minor component of ultra high energy photons. This 
component may receive contributions from different 
sources. Besides the expected flux generated by the prop- 
agation of the EECR in the intergalactic medium (T), 
photons may also be originated in different astrophysical 
environments as by-products of particle acceleration in 



nearby cosmic ray sources (e.g. [2|) and, fundamentally, 
in top-down scenarios involving the decay of super 
heavy relic particles or topological defects [3|. Extreme 
energy photons have not been unambiguously observed 
yet. However, it is expected that JEM-EUSO 0, 0, 
with its unprecedented exposure will change this situa- 
tion in the next few years. 

In the present work we estimate the photon-proton 
discrimination power of JEM-EUSO. In particular we 
develop two complementary techniques to evaluate an 
upper limit on the fraction of photons relative to pro- 
ton primaries in the integral cosmic rays flux by us- 
ing the atmospheric depth of maximum development, 
X max , of the corresponding atmospheric showers. The 
longitudinal evolution of high energy photon showers 
is dominated by the interplay between magnetospheric 
photon splitting and the LPM effect. The first process 
is highly dependent on the incoming direction of the 
photon with respect to the geomagnetic field and its 
intensity. Therefore, the translation of JEM-EUSO along 
its orbit is a distinctive parameter which adds richness 
and complexity to the analysis with respect to a tradi- 
tional Earth bound observatory and must be taken into 
account. 

II. Numerical Approach 

Given a sample of N events, an ideal upper limit to 
the photon fraction may be calculated under the a priori 
assumption that actually no photon exists in the sample: 

T™ n = 1 - (1 - a) 1/N (1) 

where a is the confidence level of rejection. However, 
in practice, the probability of the existence of photons 
must be realistically assessed through some observa- 
tional technique which involves the determination of 
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Fig. 1. Contour plot of X ma 
and longitude in the Earth. The showers considered are such that 
[30°, 60°] and E £ [10 19s \ 10 20 ] eV. 



experimental parameters, which leads unavoidably to 
less restrictive upper limits than the previous one. 

In this work the X max parameter is considered for 
the discrimination of protons and photons showers. A 
shower library was generated by using the program 
CONEX H which consist of 3 x 10 5 proton showers 
following a power law energy spectrum of spectral 
index 7 = -2.7 in the interval [10 19 8 , 10 21 ] eV 
and with uniformly distributed arrival directions. Also 
> 5.1 x 10 5 photon showers were generated under 
the same conditions but in this case cores were also 
uniformly distributed on the surface of the Earth in order 
to properly take into account pre-showering (i.e., photon 
splitting) in the geomagnetic field. 

A measure of the discrimination power of the X max 
parameter is given by the merit factor, 



V = 



(2) 



where med[X^ ax ] is the median of X max (A = j,pr) 
distribution and AX^ ax is one half of the length of the 
region of 68% of probability of X^ ax distribution. 

Fig. Q] shows a contour plot of X max — Xf y l ax as a 
function of latitude and longitude of the core on the 
Earth, for 9 € [30°, 60°] and E € [10 19 8 , 10 20 ] eV. 
It can be seen that there are regions over the Earth 
surface where this difference is larger and in which 
the discrimination between protons and photons is more 
efficient. 

Motivated by this result the concept of mask, 
Q(r)Li m ), is introduced as those regions over the Earth 
surface where 77 is larger than a given value r\Lim- 
Fig. [2] shows the median of X max and the region of 
68% of probability as a function of primary energy for 
protons and photons with 9 E [30°, 60°]. For the case of 
photons different masks are considered, Cl(r]Li m = 0) 
(all events), fl(r]Lim — 1) and fl(r]Lim — 1.5). Note 
that the masks are functions of primary energy. From 
the figure it can be seen that for the photons there are 
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Fig. 2. Median and region of 68% of probability of X 

max as a 

function of primary energy for 9 £ [30°, 60°]. In the case of photons 



different masks are considered, r;^i„ 
and rj Lim = 1.5. 
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Fig. 3. Fraction of events as a function of ?7£; m for different cuts 
in zenith angle and for E G [10 19 8 , 10 20 ] eV. 



two well defined regions, the first corresponds to primary 
energies smaller than 10 20 1 eV in which the X^ ax 
distribution is composed by LPM dominated showers 
that both, suffer and do not suffer photon splitting in the 
Earth magnetic field. In the second region all showers 
undergo photon splitting which generates, on average, 
smaller values of X^ ax and smaller fluctuations Q- 

The X? nax distributions obtained for masks with 
larger rjLi m allow a better separation between protons 
and photons. However, the total number of events also 
depends on the assumed mask and, in particular, de- 
creases with rjLim. Fig. |3] shows the fraction of events 
as a function of r\u m for different cuts in zenith angle 
and for E £ [10 19 8 , 10 20 ] eV. It can be seen that for 
vertical showers larger values of 77 are obtained which 
are also distributed over a wider range. 

Two methods were developed in order to calculate an 
upper limit for the photon fraction. The first one is based 
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on the abundance estimator first introduced in 1 8 1 



N 



/7PQ 



(3) 



where f 7 (X max ) and f pr {X max ) are the photon and 
proton distribution functions, X % max are experimental 
values of X max and iV is the sample size. £,x max is an 
estimator of the photon abundance, c 7 = N 1 /N where 
N 1 is the number of photons in the sample. The mean 
value of £,x mal is a linear function of c 7 , its standard 



deviation is proportional to 1/ v iV and, for large values 
of N, it follows a Gaussian distribution. 

For the case in which £,x max is compatible with a pure 
proton sample, an upper limit to the photon fraction, 
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The distribution functions needed to calculate c 
are obtained from the simulated data by using the non 
parametric method of kernel superposition with adaptive 
bandwidth (5). 

It can be seen from Eq. (0]i that, the larger the sample 
size, the smaller c™ m . Although, it is not obvious from 
this expression, it is possible to show that for larger 
values of 77 also smaller values of clj™ are obtained. 
Fig. |4] shows c'"" n as a function of r\Lim, i-e- for 
different masks, for E € [10 19 8 , 10 20 ] eV, 9 € [30°, 60°] 
and 9 £ [45°, 90°] and with and without assuming a 
Gaussian uncertainty of 70 g cm -2 for 9 £ [30°, 60°] 
and 60 g cm~ 2 for 6 £ [45°, 90°]. It can be seen that 
c mm i ncreases with r]Lim, which means that although 
the discriminator power of X max increases the number 
of events decreases so rapidly producing larger values 
of c™", i.e. in this case the number of events is more 
important than the discrimination power of X max for a 
given mask. 

The second method developed here consists in finding 
a cut on X max which minimizes the expression of the 
upper limit obtained assuming a pure proton composi- 
tion, 

N a (NF pr (X c max )) 



FuLyXmax) 



(9) 



NF-y {X^ ax ) 

where N a (n) is the upper limit on the number of 
photons n at a confidence level a obtained assuming 
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Fig. 4. c™ in as a function of r] Llm for E G [10 19 8 , 10 20 ] eV, 
6 G [30°, 60°] and 9 G [45°, 90°] and with and without assuming 
a Gaussian uncertainty of 70 g cm -2 for 6 G [30°, 60°] and 60 g 
cm~ 2 for 6 G [45°, 90°]. 
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Fig. 5. Estimates of the protons and photons distribution functions of 
Xmax for threshold energies from 10 19 8 eV to 10 20 3 eV in steps 
of log(_B/eV) = 0.1 and for G [30°, 60°]. A power law energy 
spectrum of spectral index 7 = —4.2 is assumed. The mean value of 
the distributions increases with primary energy. 



a Poisson distribution and 



-^M-^maa:) = / dX max fA(X max ). (10) 

In order to study the upper limit for a given threshold 
energy the distribution functions of X max for protons 
and photons are obtained from MC data, by using the 
non-parametric method of kernel superposition men- 
tioned above. Fig. [5] shows the estimates of the proton 
and photon distribution functions for 9 £ [30°, 60°] 
and for threshold energies from 10 19 8 eV to 10 20 3 
eV in steps of log(£J/eV) = 0.1. A power law energy 
spectrum of spectral index 7 = —4.2 following the 
shape of the Auger spectrum iflOl is assumed. The 
figure shows, as expected, that as the threshold energy 
increases the bump due to the photons that do not 
suffer photon splitting in the geomagnetic field becomes 
progressively less important. 

Fig. |6] shows Tuh{X^ lax ) as a function of X^ ax 
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Fig. 6. Upper limit as a function of X{L 0!C for threshold energies 
from 10 19 ' 8 eV to 10 20 ' 3 eV in steps of log(S/eV) = 0.1 and for 
6 [30°, 60°]. 

obtained by using the distribution functions of Fig. [5] for 
a = 0.95 and for a total number of events above 10 19 8 
eV N — 2250 (4500 events in total but half of them have 
9 e [30°, 60°]). It can be seen that T VL {X c max ) reaches 
a minimum which depends on the threshold energy. Note 
that there is a transition at E t h = 10 20 eV from which 
the minimum is reached at larger values of X max . This is 
due to the change in the shape of the photon distribution 
function when the threshold energy increases (see Fig. 
|5j. Finally, the upper limit is obtained evaluating J-jjl 
in X^ ax of the minimum. 

Fig. [7] shows the upper limits on the fraction of pho- 
tons in the integral cosmic ray flux, at 95% of confidence 
level, obtained in the ideal case jF" lm (dashed line), by 
using the £,x maT method for 8 € [30° , 60°] and assuming 
a 70 g cm -2 of Gaussian uncertainty and no uncertainty 
in the determination of X max (solid and dash-three dots- 
dash gray lines, respectively), by using the optimized 
cut method for 6 € [30°, 60°] and assuming a and 70 
g cm -2 of Gaussian uncertainty in the determination of 
X max (dash-dot-dash and solid black lines, respectively) 
and the upper limits obtained by different experiments. 
It can be seen that the upper limits obtained by using 
both methods introduced here are about one order of 
magnitude larger than the ideal case. This is due to the 
limitation imposed by the X max parameter to discrim- 
inate between photons and protons. Moreover, for the 
energies > 10 20 eV, the X max distribution of photon 
showers is dominated by photon splitting decreasing the 
discrimination power of this parameter. 

Fig. [7] also shows that for energies bellow 10 20 1 
eV the optimized cut method results better than £x maa 
method. This happens because the former takes advan- 
tage of the part of the distribution function originated by 
photons that do not suffer photon splitting. For larger 
energies the £x max method results better because the 
part of the X max distribution of unconverted photons is 
less important and this methods takes into account the 
whole shape of proton and photon distribution functions. 




Fig. 7. The upper limits on the fraction of photons in the integral 
cosmic ray flux at 95% of confidence level. Dashed line corresponds 
to the ideal case in which it is known that there is no photons in 
the data. Solid and dash-three dots-dash gray lines are the upper 
limits obtained by using £.x max method assuming a 70 g cm -2 
of Gaussian uncertainty and no uncertainty in the determination of 
X ma x, respectively. Solid and dash-dot-dash black lines are the upper 
limits obtained by using the optimized cut method assuming a 70 g 
cm~ 2 of Gaussian uncertainty and no uncertainty in the determination 
of Xmax, respectively. Shadow region is the prediction for the GZK 
photons [T]. Black arrows are experimental limits, HP: Haverah Park 
EH; Ai, A 2 : AGASA El, (B); A™, A SD : Auger O, 03); AY: 
AGAS A- Yakutsk Q6); Y: Yakutsk Q7J. 



III. Conclusions 

In the present work we demonstrate that, based on 
pure statistics, JEM-EUSO will be able to set an upper 
limit to the photon fraction at the highest energies of 
~ 10 _1 %, well inside the GZK-photon flux expectation. 
A more realistic estimate, taking into account a con- 
servative discrimination parameter and its experimental 
uncertainty the upper limit is larger but still considerably 
smaller than the corresponding values set so far by 
existing experiments. Therefore, given its large exposure, 
JEM-EUSO should be the first cosmic ray experiment to 
unambiguously detect cosmogenic photons. 
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